J/ijj and i] c in the nuclear medium: QCD sum rule approach* 
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Here q = (to,q), and \) n . m . is the ground state of nuclear 
matter which we take to be at rest. For the J /if we take 
the vector current jY = cry^c and for the rj c , we use the 
pseudoscalar current j = ic^y^c. In the region of large 
and positive Q 2 = q 2 — uj 2 we can express the correlation 
function through an operator product expansion (short 
distance expansion) and write the left hand side of 
eq.(0) as 
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Investigating the behaviour of heavy quark systems in 
a nuclear medium is of great interest, for several rea- 
sons. First, the ongoing discussion of J/ip suppression in 
ultrarelativisitic heavy-ion collisions as a possible quark- 
gluon plasma signal requires detailed knowledge about 
the in-medium interactions of the J/if under "normal" , 
non-plasma conditions. Furthermore, as Brodsky et al Jl]] 
pointed out, multigluon exchange can lead to an attrac- 
tive potential between a cc-meson and a nucleon, such 
that, for example, the r\ c could form bound states even 
with light nuclei. In more recent calculations the esti- 
mated charmonium binding energy in nuclear systems 
was found to be of the order of 10 MeV 

In the present paper we study the in-medium be- 
haviour of the J I if) and r\ c using QCD sum rules ||. The 
QCD sum rule approach connects the spectral density of 
a given current correlation function via a dispersion rela- 
tion with the QCD operator product expansion (OPE). 
In-medium QCD sum rules have so far been applied only 
for light quark systems, in order to study possible shifts 
of the in-medium masses of nucleons and vector 

mesons Such calculations suffer from uncertainties, 
e.g. due to assumptions about factorization of four-quark 
condensates which may not be justified. As we shall see, 
in-medium QCD sum rules applied to heavy quark sys- 
tems are expected to be more reliable. Up to dimen- 
sion four, the order to which the vacuum sum rules for 
hadrons involving heavy quarks are commonly expanded, 
all condensate parameters are quite well known and there 
are no ambiguities in the OPE. We also find that un- 
certainties caused by possibly large hadronic in-medium 
decay widths are much smaller than for light-quark sys- 
tems. 

Our starting point is the time ordered current-current 
correlation function of two heavy quark currents in nu- 
clear matter, 
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Here the O n are operators of (mass) dimension n, renor- 
malized at a scale [i 2 , and C n are the perturbative Wilson 
coefficients. 

At baryon densities Pn for which the chemical poten- 
tial is small compared to the scale fi separating short and 
long distance phenomena, all density effects can be put 
into the pjq dependence of the condensates (Oan); an d we 
can use the perturbative Wilson coefficients calculated 
in the vacuum |l0|]l2j ]. In heavy quark systems the ex- 
pansion of quark operators in terms of inverse powers of 
the large quark mass permits to express them entirely in 
terms of gluonic operators |^|,|l3],|l4| . In the vacuum only 
the scalar gluon condensate (^f-G ^G^ v ) contributes up 
to dimension four. In nuclear matter, an additional con- 
tribution involving in-medium expectation values of the 
twist-2 tensorial gluon operator (^-G aa Gp) enters. We 
discuss this new term in some detail. 

We will use the linear, low-density approximation |lq ] 
for the in-medium condensates: 
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where ()o represents the vacuum expectation value, and 
the nucleon state (taken at rest in eq. (|^)) is normalized 
as (N(p')\N(p)) = 2p (2Tr) 3 S 3 (p-p'). The in-medium 
changes of the condensates can then be related to the nu- 
cleon expectation values of the corresponding operators. 
For the traceless and symmetric gluonic twist- 2 tensor 
operator we write 

(N(p)\^G a °G^\N(p)) = ~(p a p? ~ I fl «V)— ^3 

7T 4 7T 

(4) 

where mjv is the nucleon mass and Aq is related to the 
following moment of the gluon distribution function G: 

A G (p 2 ) = 2 [ dxxG(x,p 2 ). (5) 



It represents twice the momentum fraction carried by glu- 
ons in the nucleon. We take Ac(&m 2 ) ~ 0. 9 Jlq] at the 
scale (U used previously by Reinders et al. |13| , |14| . The 
scalar gluon condensate (^-G 2 ) changes with density ac- 
cording to 



(—G^G^)n, n . = ( — G^Cno ~ ^m° N p N , (6) 



where rrfi N ~ 750 MeV is the nucleon mass in the chiral 
limit @. 

For the J/tp current, using the background field tech- 
nique [ jl8| , we find that the additional contribution aris- 
ing from the twist- 2 operator looks as follows, 
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where Jn = J dx[l + x(l — x)Q 2 /m 2 ]^ N . In the present 
work we study the cc-system at rest (relative to the sur- 
rounding nuclear matter) and set q = 0, so that eqs. (|l|j^) 
refer to the Euclidean region uj 2 — — Q 2 < 0. Then there 
is only one invariant function, 
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which reduces to the usual vacuum polarization function 
when the nuclear density goes to zero. 

Similarly, for the pseudoscalar case, the gluonic twist-2 
correction in the OPE has the following form: 
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Here we introduce the (dimensionless) polarization func- 
tion 

lF>,q = 0) 



IF"(-Q 2 =^) = 



(10) 



which reduces in the limit pjq — > to the usual vacuum 
polarization function. 

Our analysis is based on the moments of the polariza- 
tion function II J with J — V, P referring to the vector or 
pseudoscalar channels. The ra-th moment is connected, 
on the other side, with a dispersion integral involving 
Imf[ J , 
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Imn J (s) 
W (s + Qo)" +1 



ds, 



(11) 



at a fixed Q\ = Am 2 !;. Direct evaluation of these mo- 



ments using the OPE gives 
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The common factor A J n results from the bare loop dia- 
gram. The coefficient takes into account perturba- 
tive radiative corrections, while b J n is associated with the 
gluon condensate term, 
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The coefficients A^, a ,J n and b J n are listed in ref. |f]. The 
new contribution from the twist-2 gluon operator involves 



2tt 2 ^-A g 



3 (Am 2 .) 



jr^mNPN- 



(14) 



For the additional Wilson coefficient c„ we find in the 
vector channel: 



£(0 = &n(0- 
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3(2n + 5)(l + £) 2 F(n,i;n+f ;T ^) 

(15) 



and in the pseudoscalar channel we obtain 
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channel and c = 3/8-7T 2 in the pseudoscalar channel. In 
the vector channel the couplings /j of the crf^c current 
to the J/tp, ip', "0"; ■■■ resonances are determined by their 
(16) measured decay widths into e + e~ . Inserting eq. (ffq) into 



with hypergeometric functions F (a, 6; c; z). By compari- 
son with ref. (l3| we see that the c„'s differ very little from 



the 6 n 's. From the resulting term 6„( 



in eq. 11 



one then observes that the gluon condensate effectively 
changes by the following density dependent correction: 



( — G ) -> ( — G )o - -m u N + -m N — A G p N 
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G 2 )o(l-0.06/Wp ), 



(17) 



using (^G 2 ) = (0.35GeV) 4 and p = 0.17 fm~ 3 . 

The spectral function under the integral on the r.h.s. 
of eq.(p"lj) is parameterized as 

imn( s ) = fA s - m i) + cd ( s - s o) (i + ^ 

i \ " 

(18) 

in terms of a sum over low-lying resonances and a contin- 
uum part starting from sq, with c = 1/4-7T 2 in the vector 



eq. (p_lj), it is convenient to write 
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where m is the mass of the lowest state, the one of in- 
terest. The contributions of higher resonances as well 
as the continuum are absorbed in S^. Clearly, the rel- 
ative importance of these higher energy parts of the 
spectrum decreases with increasing n. It is common 
practice to take the ratio of two neighboring moments, 
M n _i/M„ = (m 2 + Q 2 )(l + <y n _i)/(l + S n ), so that /„ 
drops out and one can focus on the mass m. For n > 5 it 
turns out that (l+<5„_i)/ (l+6 n ) is close to one. Then the 
moment ratio does not depend on details of the higher 
resonances and continuum parts of the spectrum, and we 
have 
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FIG. 1. The r\ c and J/ip masses calculated according to eq. ( po| ) for different n at ^ = 1. We show the result in medium at 
pN = 0.17 fm (solid line) in comparison with the vacuum result (dashed line). 



The actual mass determination is done using moments 
in the range 7 < n < 11 and choosing £ = 1, just as in the 
vacuum case studied previously jl3|,[l4| . This range min- 
imizes the sensitivity to details of the high-energy spec- 
trum. Going to larger n would not be justified without 
introducing additional, unknown condensates of higher 
dimension in the OPE. 

In Fig. 1 we show the results for in-medium masses 
(solid lines) of the J/if> and rj c at normal nuclear matter 
density (pn = Po = 0.17 fm~ 3 ) in comparison with their 



vacuum values (dashed lines). Using a s (8m 2 ) = 0.21, 
m c = 1.24 GeV, <fib — 1-8 • 10~ 3 in the vacuum and 
(\>b = 1.7 • 10~ 3 , <f> c = -1.25 • 10~ 5 in nuclear matter, 
we find the following mass shifts taken at the minimal 
values of eq. (p0|): 



-7MeV, 
-5MeV. 



(21) 
(22) 



These shifts depend only very weakly on our choice of 
parameters. Taking higher resonances explicitly into ac- 
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count would lead to somewhat smaller mass shifts. 

The sensitivity to possible enlarged in-medium widths 
of the J/ip or r\ c turns out to be marginal. In fact no in- 
elastic channels exist for ground state cc mesons at rest 
interacting with a nuclcon. Even if such channels were 
present, they would not affect the mass shift analysis un- 
less the corresponding widths would reach magnitudes 
of 100 MeV or larger. This is in qualitative contrast to 
light quark systems, such as the p meson, for which the 
in-medium width becomes so large that the QCD sum 
rule analysis of a possible mass shift becomes ambiguous 
and inconclusive [H9p(|. 



In summary our in-medium QCD sum rule analysis, 
with the operator product expansion calculated up to 
dimension four, predicts attractive mass shifts of about 
5-10 MeV for J/tp and rj c in nuclear matter. This corre- 
sponds to small J/ip- and 77 c -nucleon scattering lengths 
a = — [i r Am/(2irpN) ~ (0.1 — 0.2) fm (p r is the meson- 
nucleon reduced mass). Our results for the mass shifts 
of the lowest cc states are surprisingly close to those re- 
ported in ref. Most of the calculated mass shift 
comes from the density dependence of the gluon conden- 
sate. The new term related to the fraction of momentum 
carried by gluons in the nucleon contributes less than 10% 
to the total effect. The influence of the decay widths is 
expected to be very small, at least for J/ip and rj c at 
rest [EL Of course, for charmonium systems traversing 
nuclear matter at high energy, the scattering amplitudes 
can develop substantial imaginary parts from reactions 
with nucleons producing open charm pl[ . 

After submission of this paper a similar calculation has 
been reported in ref. [p2| 
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